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We present a GUT-model-independent calculation of hadron matrix elements for 
all dimension-six operators associated with baryon number violating processes us- 
ing lattice QCD. Our results cover all the matrix elements required to estimate 
the partial lifetimes of (proton, neutron)— >(7r, /C, r;) e+, decay modes. We 
point out the necessity of disentangling two form factors that contribute to the 
matrix elements; previous calculations did not make the separation, which led to 
an underestimate of the physical matrix elements. With a correct separation, we 
find that the matrix elements have values 3 — 5 times larger than the smallest es- 
timates employed in phenomenological analyses of the nucleon decays, which gives 
stronger constraints on GUT models. We also find that the values of the matrix 
elements are comparable with the tree-level predictions of chiral Lagrangian. 



1 Introduction 

Nucleon decay is one of the most exciting predictions of grand unified theories 
(GUTs) regardless of the existence of supersymmetry (SUSY). Although none 
of the decay modes have been detected up to now, experimental effortsjDver the 
years have pushed the lower limit on the partial lifetimes of the nucleoicl, which 
can give a constraint on (SUSY-)GUTs. Moreover, exciting plans for the next 
generation of the Super-Kamiokande £xperiment to search the nucleon decay 
events are proposed in this conferenccl. 

On the theoretical side predictions of the nucleon partial lifetimes suffer 
from various uncertainties. In general the partial lifetime of the nucleon decay 
process N PS -\- 1, where N, PS and I denote nucleon, pseudoscalar meson 
and antilepton respectively, is given by 

t^\{PS\0^\N)\-'-\Fgvt\-\ (1) 

with the baryon number violating operator that appears in the low-energy 
effective Lagrangian of (SUSY-)GUTs. Fqut denotes some function in terms 
of parameters defined in the (SUSY-)GUT models. Although the precise num- 
bers of the nucleon partial lifetimes depend on the details of the theories be- 
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yond the standard model, which contain many unknown parameters, another 
main source of uncertainty is found in the evaluation of the hadronic part 
{PS\0^\N) in eq. (0). The matrix elements have been estimated by employ- 
ing various QCD models. Their results, however, scatter over the range whose 
minimum and maximum values differ by a factor of tercl. Therefore, a precise 
determination of the nucleon decay matrix elements from the first principles 
using lattice QCD is of extreme importance. 

Most important feature of our calculation is GUT model independence. All 
dimension-six operators associated with baryon number violating processes are 
classified into four types UWier the requirement of SU(3)xSU(2)xU(l) invari- 
ance at low-energy scaleEI'B. If one specifies the decay processes of interest, 
namely the processes among (proton,neutron)^(7r, i^, 77) -\-{i>,e^ , ^'^), we can 
list a complete set of independent matrix elements in QCD, and we calculate 
all the matrix elements. 

We employed two methods to evaluate the nucleon decay matrix elements 
in lattice QCD. One is the indirect method in which the matrix element 
{PS\0^\N) is estimated from (0|C''^|iV) measured on the lattice employing 
the tree-level results of chiral Lagrangian. The other is the direct method 
which directly measures (PS'lOfJiV) on the lattice. In the previous lattice 
QCD studies the indirect methocBfl gave values of the matrix elements 2 or 3 
times smaller tlian those obtained by the direct methocEltl. Recently, however, 
we pointed ouiij that the nucleon decay matrix elements {PS\0^\N) allow the 
contributions of two form factors for general lepton momentum, which should 
be disentangled in the direct method, and this explains the discrepancy be- 
tween the direct,|and indirect estimations of the matrix elements found in the 
previous studieslu where the separation was not made. 

In this report we first formulate the method to calculate the nucleon decay 
matrix elements using the lattice QCD in Sec. ||. Our results and conclusions 
are given in Sees. || and ^. More details are found in Ref. 10. 



2 Formulation of the method 

2.1 Independent matrix elements for nucleon decays 

The low energy effective theory in the baryon number violating processes is de- 
scribed in terms of SU(3)xSU(2)xU(l) gauge symmetry based on the strong 
and the electroweak interactions, which enables us to make a GUT-model- 
independent analysis. Our interest is focused on the dimension-six operators 
which are the lowest dimensional ones in the low energy effective Hamiltonian: 
operators must contain at least three quark fields to form SU(3) color singlet, 
and then an additional lepton field is required to construct a Lorentz scalar. 
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Higher-dimensional operators are suppressed by inverse powers of heavy par- 
ticle mass characterized by the theory beyond the standard model. 

All dimension-six operators are classified Jnto the four types under the 
requirement of SU(3)xSU(2)xU(l) invariance qB: 

^ilL = i^iaRUjbR)iQakcLL/3dL)£ijk<^af3, (2) 

(2) ~c ~c 

^abcd = {QaiaLQpjbL){UkcRLdRyijk^aP, (3) 

(3) ~c ~c 

^abcd = {QaiaLQ0jbL){Q'ykcLL5dL)(^ijk<^aS<^l3'y, (4) 

^ItL = 0taRUjbR)iUkcRLdR)f^ijk, (5) 

where ^p'^ = tp'^C with C = 7472 the charge conjugation matrix; i, j and k are 
SU(3) color indices; a, (3, 7 and S are SU(2) indices; a, b, c and d are generation 
indices; Ll and Ql are generic lepton and quark SU(2) doublets with the left- 
handed projection Pl = (1 — 75)/2; La, Ur, and Da are generic charged lepton 
and quark SU(2) singlets with the right-handed projection Pr = (1 +75)72. 
Fierz transformations are used to eliminate all the vector and tensor Dirac 
structures in eqs. (^ — (||). 

Our interest exists in the decay processes from the nucleon to one pseu- 
doscalar meson: (proton, neutron)^{7r, iiT, 77) -!-(/>, e"*", ^"'"). Once these decay 
modes are specified, we can list the set of independent matrix elements in QCD 
from the operators of eqs. (H) — (||): 

{TT''\e,Mu'^CPR,Ld')PLu''\p), (6) 

{7T+\e,,k{u'^CPR,Ld')PLd''\p), (7) 

{K°\e,,k {u'^CPb„ls')Plu'' \p) , (8) 

{K+\e,,k {u'^CPR^LsnPLd" \p), (9) 

{K+\e,,k{u'^CPR,Ld^)PLs''\p), (10) 

{K''\e,,k{u'^CPR,Ls')PLd''\n), (11) 

{vh,k{u'^CPR,Ld^)PLu''\p), (12) 

where we assume SU(2) isospin symmetry m„ = md and use the relations 

{PS\Olr.\N) = {PS\Orl\N), (13) 
{PS\Orr\N) = {PS\Oll\N), (14) 

due to the parity invariance. All we have to calculate in lattice QCD are 
these 14 matrix elements. Other matrix elements can be obtained by using the 
exchange of the up and down quarks and the relations of eqs. ([l^ ) and (p^. 
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2.2 Form factors in nucleon decay matrix elements 

Under the requirement of Lorentz and parity invariance, the matrix elements 
between the nucleon and the pseudoscalar meson in eqs. (^) — (|l^) can have 
two form factors: 

{PS{p)\Of\N^^\k)) = Pl {Wo{q') - W,{q')ti) u^-^\ (15) 

where Oj^ represents the three-quark operator projected to the left-handed 
chiral state, u*-^-* denotes the Dirac spinor for nucleon with either the up (s = 1) 
or down (s — 2) spin state, and is the momentum squared of the out- 



going antilepton. The contribution of the Wg term in eq.(15) is negligible 
in the physical decay amplitude, because its contribution is of the order of 
the lepton mass m; after the multiplication with antilepton spinor. However, 
since the relative magnitude of the two form factors Wq and Wq is a priori 
not known, we have to disentangle these two form factors in the lattice QCD 
calculation. Hereafter we refer to Wq and Wq as relevant and irrelevant form 
factor respectively. 

In our lattice QCD calculation we choose /c = for the nucleon spatial 
momentum and p = k — q ^ for the PS meson. In this case the Dirac 



structure of the right hand side in cq. (15) is given by 

(Wo - Wql^) U(^) 

^-!'^^^%,7^«.^^\l').(^) (16) 
Wqq -a Wo+ iqiWq J 



Wo + (mjv - x/m''ps+P^)Wq Wqp- 



-Wqp- a Wo - [niN - y/m%g+p^)Wq 

where Wo — Wqi^ is expressed by a 2x 2 block notation; a are the Pauli matrices, 

and = (1,0,0,0) or (0, 1,0,0). Using the (1, 1) and (2, 1) components in 

the 2x2 block notation of eq. (|l^), where the other components vanish, we 
can extract the relevant form factor Wo- 

The need for the separation of the contribution of the irrelevant fq 
factor was not recognized in the previous studies with the direct methoc 
The values found in these studies correspond to Wo — iqiWq instead of Wo- 
We examine how much this affects the estimate of the matrix elements in 
Sec. I 

2-3 Calculational methods 

The nucleon decay matrix elements of eq. ( |l5| ) are calculated with two methods 
referred to as the direct and the indirect ones. The former is to extract the 
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matrix elements from the three-point function of the nucleon, the PS meson 

and the baryon mimbcr violating operator. The latter is to estimate them with 
the aid of chiral Lagrangian, where we have two unknown parameters to be 
determined by the lattice QCD calculation. 

In the direct method we calculate the following ratio of the hadron three- 
point function divided by the propagators of the pseudoscalar meson and the 
nucleon: 

R{t, t') 



^'^-^^'-^^ {Jpsix', t')4si^, *)) Ex(-/iv,«(x, i) Jjv,«(0)) 

-r^{PS{p)\dfjN^^\k = 0)) t' » t » 0. (17) 

Here (5^ ^ denotes the rcnormalized operator in the naive dimensional reg- 
ularization(NDR) with the MS subtraction scheme, and Jps and Jjv,s are 
interpolating fields for the PS meson and the nucleon, respectively. The sub- 
scripts 7 and s are spinor indices; we can specify the spin state of the initial 
nucleon at rest by choosing s = 1 or 2. Zps and represent the residues of 
the PS meson propagator and the nucleon propagator. L^LyL^ is the spatial 
volume of the lattice. In lattice QCD calculation we evaluate the Green func- 
tions numerically by the Monte Carlo method with supercomputer. We move 
the baryon number violating operator ^ in terms of t between the nucleon 
source placed at t = and the PS meson sink fixed at some t' well separated 
from i = 0. Under the condition that the baryon number violating operator 
is sufficiently distanced from both the PS meson field and the nucleon field, 
which is required to avoid excited state contaminations, we can extract the 
matrix elements normalized by the spatial volume of the lattice. 

The indirect method uses the tree-level results of chiral Lagrangian, which 
are given by 



{A{udR)UL\p) 


= aPLUp[l + D + F\/{s/2f), 


(18) 


{A{udL)UL\p) 


= l3PLUp[l + D + F]/{s/2f), 


(19) 


{T:+\{udR)dL\p) 


= aPLUp [1+D + F] If, 


(20) 


{T,+ \{udL)dL\p) 


= pPLUp + D + //, 


(21) 


{K'>\{usr)ul\p) 


= aPLUp [-1 - (D - F)mN/B] //, 


(22) 


{K\usl)ul\p) 


= (iPLUp [l-iD- F)mN/B] If, 


(23) 


{K+\{uSR)dL\p) 


= aPLUp [2DmN/B] /(3/), 


(24) 
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{K+\{uSL)dL\p) 


= (3PlUp 


[2DmN/B 


] /(3/), 


(25) 


{K+\{udn)sL\p) 


= aPLUp 


[3 + 


{D + 


3F)m^./B] /(3/), 


(26) 


{K+\{udL)sL\p) 


= PPhUp 


[3 + 


iD + 


3F)mN/s] /(3/), 


(27) 


{K''\{usR)dL\n) 


= aPLUn 


[-3 


-{D 


-3F)m^/B] /(3/), 


(28) 


{K\usL)dL\n) 


= (iPLUn 


[3- 


(D- 


iF)mN/B] /(3/), 


(29) 


{r]\{udR)uL\p) 


= aPLUp 


[-1 


-{D 


- 3F)] /(^^/), 


(30) 


{r]\{udL)uL\p) 


= PPLUp 


[3- 


(D- 


3i^)]/(x/6/), 


(31) 



where / is the pion decay constant; F and D parameters are determined from 
experimental results of the semileptonic baryon decays; m^si/B = '^nI'^b 
with niB = ~ mA. We use {PS\{i\)\'^2R^L)^-iL\^) a shortened form of 
(P5|eyfe(^f CPfl„LV'^)FLVt|7V). 

The expressions in eqs. (p^— (|3l|) contain two unknown coefficients a and 
/3, which are to be determined by lattice QCD calculation. The definitions of 
a and/3 parameters are given by 

iS)\e,,k{ifcPRd^)PLU%^'^)^ciPLvS'\ (32) 
(0|e,,fc(u^^CPid^)Fi^i'=|p(^)) =/3Piw(^), (33) 

where operators are renormalized in the NDR scheme, p denotes the proton 
state. These matrix elements are obtained from the ratio of two-point func- 
tions: 

^a/3.,^ ^ Y.sk^^,u{^^CPR^Ld^)PLu\x,t)Jp,M) 

^' j:,{JpA^,t)jpAo)) ^ 

^ {0\e,,kiu'CPR^Ld')PLu''\p^'^) t»0. (34) 

Incorporating the a and (3 values determined by the lattice calculation into the 
tree-level results of chiral Lagrangian, we can obtain the values of the nuclcon 
decay matrix elements in eqs. (0) — (12). 



3 Results 

3.1 Details of numerical simulation 

Our calculation is carried out with the Wilson quark action in quenched QCD 
at /3 = 6.0 on a 28^ x 48 x 80 lattice. Gauge configurations are generated 
with the single plaquette action separated by 2000 pseudo heat-bath sweeps. 
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R(t,t'=29) for <it l(ud|,)Uj^lp> 



0.00 


1 1 1 1 1 1 1 
<7t°l(udg)U|^lp> 






-0.02 
-0.04 


i 



0.O4 0.08 

2 2 

-q a 



Figure 1: Ratio B.{t,t' = 29) for the rel- 
evant form factor Wq in {iT^\{udii)ui^\p) . 
Solid lines denote the fitting results with 
an error band of one standard deviation. 



Figure 2: —q a dependences for Wo in 
{TT^\(udfi)ui;^\p). Combination of form fac- 
tors Wo — iq^Wq is also plotted for com- 
parison. Solid line denotes the function 
CO + ci ■ {-q^a?) + C2- 



-q^a'^f. 



We analyzed 100 configurations for the calculation of the nucleon decay ma- 
trix elements after the thermalization of 22000 sweeps. The four hopping 
parameters K = 0.15620, 0.15568, 0.15516 and 0.15464 are adopted such that 
the physical point for the K meson can be interpolated. The critical hop- 
ping parameter Kc — 0.15714(1) is determined by extrapolating the results 
of at the four hopping parameters linearly in 1/2K to = 0. The p 
meson mass at the chiral limit is used to determine the inverse lattice spac- 
ing a^^ — 2.30(4)GeV with mp — 770MeV as input. The strange quark mass 
iTLsa = 0.0464(16)(ifs — 0.15488(7)), which is estimated from the experimental 
ratio TTiK/mp = 0.644, is in the middle of = 0.15516 and K = 0.15464. 

We calculate the ratio of eq. (17) with the nucleon field fixed at t = and 
the PS meson field aXt = 29. Four spatial momenta pa = (0, 0, 0), (7r/14, 0, 0), 
(0,7r/14,0) and (0,0, 7r/24) are imposed on the PS meson in the final state. 
For the p 7^ cases we provide different quark masses for the valence quark 
connecting Jps and ^ and the other valence quarks: TO2 is for the former 
and mi for the latter. By this assignment we can distinguish the strange quark 
mass from the up and down quark mass. 



3.2 Nucleon decay matrix elements with direct and indirect methods 

We first present the results of the nucleon decay matrix elements obtained by 
the direct method. In Fig. |l| we show time dependence of R{t, t' — 29) with 

|p|a — 7r/14 for the matrix element {'K'^\eijk(u'^'^ CPiid^)PLU^\p) as a represen- 
tative case. The result of constant fit is depicted by the set of three horizontal 
lines. We choose the fitting range to be 8 < < < 16 for all the matrix elements 
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of eqs. dq)— (p^ such that the excited state contaminations in the nucleon and 
PS meson states can be avoided simuhaneously. 

Figure Ushows —q^a? dependence of the relevant form factor Wo(g^) in the 
matrix element (7r°|eijfc(u* CPRd^)PLu'^\p) ^ where the operators are renormal- 
ized with the NDR scheme at fi = 1/a. The fitting result in Fig. |l| corresponds 
to the data at -q^a^ = 0.0259(22) in Fig. |. The combination Wq - iq^Wq 
is also plotted in Fig. || for comparison, which is obtained by following the 
method in Ref. 8. The magnitude of Wo((?^) is more than two times larger 
than that of Wo{q^) - iq^Wqiq"^). 

To interpolate the relevant form factor Wq to q^ — 0, which is the on-shell 
point of outgoing antilepton, we employ the following fitting function: 

Co + ci • (-(?^) + C2 • (-(?^)^ + C3 • TOi + C4 • m2, (35) 

where we assume that the form factor could have the mi and TO2 dependences 
through the nucleon and PS meson masses. We extrapolate mi and 1112 to 
the chiral limit for the matrix elements of eqs. (§), (|) and (|l|), while 7712 IS 
interpolated to the physical strange quark mass with mi taken to the chiral 
limit for the matrix elements of eqs. (H)— (pT|). The solid line and the open 
circle at ~q^a^ = in Fig. || denote the fitting result of the data employing 
the function of eq. (Bq), where we find that the charged Icpton masses m^a^ — 
4.9 X 10~* and m^a^ = 2.1 x 10~^ are negligible in the current numerical 
statistics. The solid line expresses the function cq + ci • (— g^) + C2 • {—q^Y 
employing the fitting results of cq, ci and C2. 

Let us turn to the results obtained by the indirect method, which uses 
the tree- level results of chiral Lagrangian in eqs. (p^) — (pl|) . The a and /? 
parameters at each hopping parameter are extracted from a constant fit of 
the ratio of eq. (|3^). Applying linear fits to the data as a function of quark 
mass 7nqa = [l/K - l/Kc)/2, we obtain q;(NDR, 1/a) = -0.015(l)GeV3 
and /?(NDR, 1/a) = 0.014(l)GeV^ in the chiral limit with the use of a"^ = 
2.30(4)GeV. 

In phcnomenological QlIT model analyses of the nucleon decays, the val- 
ues |q;| = |/3| = O.OOSGeV^ liil are conservatively taken as these are the smallest 
estimate among various QCD model calculationfl. The previous lattice calcu- 
lations, however, indicated values of thfise parameters considerably larger than 
the minimum model estimate abovfllUfl. Our results, significantly improved 
over the previous ones due to the use of higher statistics, larger spatial size, 
lighter quark masses and smaller lattice spacing, have confirmed this trend: 
the values we obtained are about five times larger than \a\ = \(3\ — O.OOSGeV^. 

In Fig. ^ we compare the nucleon decay matrix elements obtained by the 
direct method with those by the indirect one using the tree-level results of 
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-<jr''l(ud^)uLlp> 
<Jt''l(udL)Uj^lp> 
-<Ji^l(udjj)dj^lp> 
<Jl*l(udL)dj^lp> 
<K°I(us^)UlIp> 
<K°l(us^^)u^^lp> 
-<K*l(us^)dj^lp> 
<K*l(usL)dLlp> 
-<K*l(udg)sLlp> 
<K*l(udj^)sLlp> 
<K°l(us^)dLln> 
<K°l(usL)dLln> 
<Til(ud^)uJp> 
<ril(udJuLlp> 

0.0 0.1 0.2 0.3 

GeV 

Figure 3: Comparison of relevant form factors with tree-level predictions of CliPT. Crosses 
denote the ChPT results with \a\ = = O.OOSGeV"^. For numerical values see Ref. 10. 

chiral Lagrangian (squares), where we employ the expressions of eqs. 
with a(NDR, 1/a) = -0.015(l)GeV3, /3(NDR, l/a) = 0.014(l)GeVyj^ = 
O.iaiGeV, WIN = 0.94GeV, ms = 1.15GeV, D = 0.80 and F = O.4M We 
observe that the two set of resuhs are roughly comparable. This leads us to 
consider that the large discrepancy between the results of the two methods 
found in Rcfs. 8,9 is mainly due to the neglect of the Wq{q'^) term in cq. (p^. 

It is also intriguing to compare our results with the tree-level predictions 
of chiral Lagrangian with |a| = \(3\ = O.OOSGeV^ (crosses). Our results with 
the direct method are 3 — 5 times larger than the smallest estimates except 
{ri\{udii)uL\p) . Hence they are expected to give stronger constraints on the 
parameters of GUT models. 

4 Conclusions 

In this article we have reported progress in the lattice study of the nucleon 
decay matrix elements. In order to enable a GUT-model-independent analysis 
of the nucleon decay, we have extracted the form factors of all the indepen- 
dent matrix elements relevant for the (proton, neutron)^ (tt, K, rj) +(P, e+, /i+) 
decay processes without invoking chiral Lagrangian. 

We have also pointed out the necessity of separating out the contribution 



X ChPT(a ,B ) 

^ min'~niin' 

□ ChPT(ours) 
• W„(q'=0) 

I I , I 
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of an irrelevant form factor in lattice calculations for a correct estimate of the 
matrix elements at the physical point. With this separation, the matrix el- 
ements obtained from the three-point functions are roughly comparable with 
the tree-level predictions of chiral Lagrangian with the a and /3 parameters de- 
termined on the same lattice. The magnitiide of the matrix elements, however, 
are 3 to 5 times larger than those with the smallest estimate of a and /? among 
various QCD model calculations. Our results would stimulate phenomenologi- 
cal interests as the larger values of the nucleon decay matrix elements can give 
more stringent constraints on GUT models. 

Major systematic errors conceivably affecting our present results are the 
scaling violations and the quenching effects. The former can be investigated by 
repeating the simulation at several lattice spacings; the latter is eliminated once 
configurations are generated with dynamical quarks, where it is straightforward 
to apply our method. We leave these points to future studies. 
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